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Abstract 

We consider the effect of weak nonmagnetic disorder on two chains of in- 
teracting fermions (with and without spins) coupled by interchain hopping. 
For the spinless case, interchain hopping increases localization for repulsive 
interactions but stabilizes the s-wave superconducting phase for attractive 
interactions. For the case with spin, the d-wave phase arising from purely 
repulsive interactions in the clean system is destroyed by an infinitesimal dis- 
order while for attractive interactions, the s-wave superconductivity is more 
resistant to disorder than in the one-chain case. 
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One dimensional electronic systems are known to be the simplest realizations of non- 
fermi liquids, and to have generic properties known as Luttinger liquids Despite the 

good understanding of purely one dimensional systems, the effects of interchain hopping, 
allowing to go from one to higher (two or three) dimensions are much less known. Whether 
such hopping is relevant and drives the system towards a Fermi liquid fixed point or whether 
non-fermi liquid properties can be retained even in presence of finite hopping is still a contro- 
versial issue . A theoretical understanding of quasi one-dimensional strongly correlated 
electronic systems (especially the crossover from Luttinger to Fermi liquid) would be rele- 
vant for the physics of organic conductors and may perhaps give some insights for High-Tc 
superconductors. Consequently, there has been in the recent years, a growing interest in 
systems of coupled interacting electron chains, and in particular in systems of two chains. 
They present the advantage to allow a careful study of the effects of the hopping, as well as 
to be tractable by powerful analytical IHbj and numerical techniques |15|-[T8|. In addition 



there exists good experimental realizations of such systems. For example Sr n _ 1 Cu n+ i0 2n 
[TJ| and VO2P2O7 |f2U| compounds are very good realizations of coupled spin chains. Upon 



doping, such compounds will give coupled Hubbard chains. Although the complete phase 
diagram of such systems is still under study, a generic property of two coupled chains system 
is the appearance of a c?-wave like superconducting phase for repulsive interactions. 

In this work we study the effects of non-magnetic disorder on such two chains systems, 
both for the case of spinless electrons and for electrons with spins. Such a study has a double 
interest: in experimental systems, disorder will be present, and it is therefore essential to 
know the stability of the phases found in the pure system. It is now well known that for 
a strictly one dimensional system, disorder has extremely strong effects and an arbitrarily 
weak disorder destroys superconductivity except for exceedingly attractive interactions [|2T] . 



In addition, on a more theoretical level, the two chain problem is the simplest one to study 
the effects of interchain hopping onto the Anderson localization in presence of interactions, 
giving some clues to the unsatisfactorily understood physics of such transition in more than 
one dimension. We show here that for the spinless model the superconducting phase for 



attractive interactions is stable towards weak disorder at the opposite of what happens for 
a purely one dimensional system. For the model with spins and attractive interactions an 
arbitrarily weak disorder destroys the superconductivity if the interactions are not attractive 
enough as in the one chain case. Nevertheless, the threshold in interaction strength to induce 
superconductivity is much smaller for two chains. In particular it can now be reached for a 
pure Hubbard attraction at variance to the one chain case f22fl . On the other hand the c?-wave 



type superconductivity found for repulsive interactions is completely unstable with respect 
to arbitrarily weak disorder. In two spinless chains, attractive interactions reduce localization 
compared to the case of a single chain whereas repulsive ones enhance localization. For the 
case with spin, two chains are always less localized than their one chain counterpart. For 
each case we also compute physical quantities such as localization length and conductivity. 

Let us consider first two chains of spinless fermions coupled by an interchain hopping t±. 
Such model can also be mapped to two spin chains coupled by an exchange X-Y term, in 
the presence of a magnetic field. For simplicity we will just consider here a nearest neighbor 
interaction V. More complicated interactions can be considered without changing the main 
physical results. Details will be given elsewhere . The disorder is modelled by a random 



on-site potential e^ p uncorrelated from site to site and from chain to chain. The Hamiltonian 
then reads 



H=-tJ2 c i, P c i+i,p + h-c. + V n itP n i+ltP 



hp » 

+t± C l,l c i,-1 + h - c - + J2 e i,P n i,P (!) 



where p = —1,1 is the chain index and i is the site index. 

It is convenient to rewrite the Hamiltonian in a standard boson representation p|-§]]- We 
therefore linearize the fermions dispersion relation around hp, introduce right movers (R) 
and left movers (L) for each chain. Then we take the continuum limit c„ jr . iP — > y/aip r ^ p (na) 
with r = L, R, p = ±1 the chain index and a the lattice spacing. We rewrite the hamiltonian 
in the bonding ip Q = and antibonding ip n = ^'Jtf 1 bands base and introduce the 

densities p r ,o^{ x ) ='■ Vr 7r( a; ) ? /v,o,7r(2 ; ) : - We then define the canonically conjugate fields <fi p u 



and IF 11 via : 

rill 



9 x (j> p ,\ 

n. 



7T 



(2) 
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(PiJ,o ± PR,n ~ PL,o =F Pi,. 



and the field 0^ ii(#) = /foo n p u(x')dx'. In term of these fields the Hamiltonian becomes for 



the pure case (ej )P = 0) : 
H = H p + H ll , 

m = 



Z7T 



7/ 

« p A- p (7rn p ) 2 + -f(^ 
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/" dxtx—d x <P {l + [ dx -^^008(7801,) + cos^y 
7 it J [ZTiaY (lixaY 



The expressions of the K, u, g in terms of the original parameters of the hamiltonian can 
easily be obtained [T^J2^]. For the pure t-V model one has K p < 1 (resp. K p > 1) and 
gf < (resp. gf > 0) for repulsive (resp. attractive ) interactions and K\\ = 1 for all t,V. 
By adding interchain interactions, one has access to the cases K p > 1 and gf < or K p < 1 
and gf > 0. The complete phase diagram in the pure case has been obtained in |14[ by a 
mapping on a problem of one chain of fermions with spin and spin-anisotropic interactions 
in a magnetic field f24|. The t± term suppresses cos(v^80||) so that 6\\ develops a gap and 
acquires a non-zero expectation value determined by minimizing the ground state energy. 
The operators with divergent associated susceptibilities are then: 



O 



CDW 71 



^RoWlAx) + ^R^L.oix) ~ costV^i 



o 



OAF 



~ e 



sin 



O 



si 



O 



S2 



Sill 



COS 



(v^ll 



These operators describe respectively out of phase charge density waves, an orbital antifer- 
romagnetic phase and chain symmetric "s" and chain antisymmetric "d" type superconduc- 
tivity. 

For gf < we have (9«) = giving an S2 phase for K p > 1 and the CDW 71 for K p < 1. 
For gf > we have (6*11) = -j= giving the SI phase for K p > 1 and the OAF phase for 
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Kp < 1. In |Tj] the bosonized forms of Osi and O52 are exchanged due to the neglect 
of anticommuting operators, so that the two superconducting phase have been erroneously 
exchanged. 

Now, we consider the effect of the disorder. Taking the continuum limit for the on-site 
random potential, keeping only the 2kp terms in the bosonized expressions (as the forward 
scattering does not induce localization ||25||), and finally going to bonding and antibonding 
bands, one finds that the coupling to disorder is represented by two terms: 

dx 



H s = J JL^( x )e^" cos(V20||) + h.c. (4) 

H a = f —U^y^" cos(V20|,) + h.c. (5) 
J na 

where £ SjQ are two uncorrelated gaussian distributed random potentials such that 



£n(x)£n'( x ') = D n 5 n:Tl i5(x — x') with n,n' = s,a. In the original lattice problem, the role of 
£ s>a would be played respectively by e„ ; i±e n _i. We consider in the following a disorder weak 
enough not to destroy the gaps opened by the interchain coupling in the pure system. This 
corresponds to the limit D <C t±. The other limit where both the interchain hopping and 
the disorder are small compared to the other parameters by for arbitrary magnitude is only 
important in the vicinity of the noninteracting point. It can be studied by similar meth- 
ods and will be discussed elsewhere [P3[ . In the weak disorder limit, <fiu has huge quantum 
fluctuations, and consequently D s will always be less relevant than D a . We can therefore 
concentrate on the latter and forget about the former. 

First, we consider the case of gj < (i.e. V > for the t-V model). In that case, we can 
replace cos(v20|i) by its (non-zero) mean value. Then the coupling to disorder fl5|) reduces 
to C J dx^a{x)e ty ^^ p<yX ^ + h.c. and the RG equations for that problem have been derived and 
analysed in [^TJ. In particular the disorder will grow under renormalization as 

^ = D a (3 - K p ) (6) 
where / = ln(a) is the standard logarithmic scale associated with cutoff renormalization. (EH 



implies a localization-delocalization transition |y] at K p = 3. As a consequence, the d-wave 



superconducting phase is unstable in the presence of disorder except for extremely strong 

attractive interactions. For a simple t-V model for which K p < 1, the CDW ground state is 

also unstable to disorder. The localized phase is a Pinned Charge Density Wave phase, with 

i 

a localization length given by L2 C h. = (1/D) 3 ~ k p . This is to be compared to the localization 

i 

length of a one dimensional spinless system L\ c h_ = (1/ D) 3 ~ 2K p . For repulsive interactions 
the effects of the interchain hopping is therefore to decrease the localization length and to 
make the two chains system more localized. The conductivity above the pinning temperature 
u/L2ch. can be obtained by methods similar to |2l| and varies as cr(T) ~ t 2 ~ k p^ 

On the other hand if one considers gj > 0, i.e. attractive interactions for a t-V model, 
then (#|| ) = -j= and in a first approximation (cos(v / 2#||)) = so that there is apparently 
no coupling at all to the disorder. Obviously, such an approximation is too crude and we 
must take into account the fluctuations of 9» around its mean value. Keeping only the 
relevant terms, and integrating out the fluctuations of 6\\ around its mean value we obtain 
the following effective action for (f> p : 

S e /f = J dxdr 

with £ eff (x)£* ff (x') = D eff 8(x - x') and D eff ~ D\. 

The renormalization of the disorder will again be given by an equation similar to (||), 
but with a coefficient (3 — 4K p )/2 in front of D a . The disorder is now relevant only for 
K p < 3/4, leading to three different phases for gf > 0: a random orbital antiferromagnet for 
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K p < 3/4 (with a localization length L 2c h. — (1/ D) 3 ~ 4K p ), an ordered orbital antiferromagnet 
for 3/4 < Kp < 1 and a s-wave superconducting phase for K p > 1. For the t-V model, 
K p > 1, and the "s"-wave superconducting phase is therefore stable with respect to weak 
disorder, at variance to the single chain problem. For the latter the derealization only 
occured for extremely attractive interactions i.e. K p > 3/2. For the two chains problem the 
localization-delocalization transition arises in the immediate vicinity of the non-interacting 
point. Contrarily to the case of repulsive interactions, interchain hopping is now strongly 
reducing the localization effects compared to the one dimensional case. The determination 
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2wK n 



+ Z eff (x)e 



lV8(f>p(x,T) 



x e 



(7) 



of the critical properties at the boundary between the repulsive (localized) regime and the 
attractive (superconducting) one, requires to treat the case where the gaps induced by the 
hopping and the disorder have arbitrary relative strength |23[] . The conductivity now behaves 



as o-(T) ~ T 2 - 4K ", and diver ges as T — > since the ground state is superconducting. In 
addition, since the disorder is less relevant for attractive interactions than for repulsive ones 
one can also expect the charge stiffness |26] for a disordered finite length two chains system 



to be larger for the attractive case than for the repulsive one, similarly to the one chain 
27] , 22] system, but with much more dramatic effects. 



Let us consider now the problem with spins. Here again, we will for simplicity only 
consider the case of a local Hubbard interaction. More general interactions can be treated 



by the same method, and give rise to a richer phase diagram ||11|| . The Hamiltonian is now 

H = —t ^ cl + i ap Ci yC r )P + h.c. — t± ^ c\ ap Ci iCr) -p 

+U n i,i,p n i,l,P + ^ e^pTii^p (8) 

i,p i,c,P 

The pure system can again be studied by using a boson representation. One introduces 
similar fields than in ([2]) for each spin degree of freedom, and make the symmetric (charge) 
<pp = 0t + 0J. an d antisymmetric (spin) CT = 0j — 0j linear combinations. One ends 
with four bosonic fields instead of two for the spinless case. The bosonized Hamiltonian 
is quite lengthy and will not be reproduced here for reasons of space. It can be found in 
|TT| , and we will use in the following the notations of this paper. All physical quantities 



depends on a parameter K p+ of the symmetric charge mode, analogous to the K p of the 
spinless problem. For the purely repulsive case, U > 0, only one of the four bosonic fields 



that describe the low-energy physics of the system (0 P +) is gapless [IT]. The mean values 



of the 3 other fields are determined by minimizing the energy of the ground state giving 



(0 p -) = 0, (4> a +) — f > ((pa-) = f j leading to a d-wave superconductive phase [TTJ. For the 
attractive case U < 0, <p p+ is again massless in the pure case. But now, we have (0 P -) = 0, 
(4> a +) = 0, (4>a-) = 0. Here, the most divergent fluctuations are associated with the operator 
Osc s ~ e ^ p+ cos(0 (7+ ) cos(0o-_) which is the order parameter for s-wave superconductivity. 



The coupling to disorder arises again via two terms 



H a = j ia{.x)0 C DW-i.x) + C( X )°CDW^( X ) dx ( 9 ) 

H s = I Ux)Ocdw°{x) + e s (x)Oh DWo (x)dx (10) 



Where £ n {x)£n'{ x ')* — D n 5 n)n '5(x — x')(n,n' = a,s), the £ n being random gaussian dis- 
tributed potentials. The operators O represents charge density waves : CDW° is the in- 
phase charge density wave, and CDW n is the out of phase one. 

Assuming again that the disorder is weak enough not to destroy the gaps, the O operators 
have the simple form for repulsive interactions 

Ocdw ~ e l *' + sin(0 p _) (11) 
Ocdw ~ sin(^_) (12) 

These two operators have exponentially decaying correlation functions and no direct coupling 
with disorder would exist if one just took into account the mean values of the fields </> Pi _ 
and Q a -. As in the spinless case one should integrate over fluctuations to get the effective 
coupling 



odisorder 



/ U.{x)e l2 ^+^ T) dxdT + h.c. (13) 



One can also view (|13|) as the coupling of the fermions with the kp ± &Ftt Fourier com- 
ponent of the disordered potential. The problem has in fact been reduced to a problem 
of spinless fermions. The localization-delocalization would occur at K p+ = 3/2 but purely 
repulsive interaction imply K < 1. The c?-wave phase is therefore unstable to abitrar- 



ily weak disorder. The symmetric ( |T0| ) and the antisymmetric (|9D part of the disorder 
contribute equally to destroy the d-wave superconductivity, in contrast with the spinless 
case where the antisymmetric part was the most relevant. The localization length in 
that phase is L-ich. ~ (1/ ' D) 2 ^ 3 ~ 2Kp ' + \ and therefore longer than the corresponding one 
L\ch. ~ (1/ ' D) 1 '( 2 ~ k p> + > of the one dimensional spinning chain. The two chains problem is 
less localized than the corresponding one dimensional one even for repulsive interactions, in 
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contrast with the spinless case. This is in qualitative agreement with what one expects in 
the absence of interactions where the localization length is proportional to the number of 
channels in the system. 

For the attractive case, the O operators take a different simplified form, due to the 
different gaps in the system 

O CW o~e^ + cos(0 p _) (14) 
Ocdw* ~ sin(0 ff _) sm(0 CT+ ) (15) 



By substituting in (H) and (|T0|) and integrating over fluctuations we end with an action of 
the form (|13|). This time, K p>+ > 1 , so that we can attain the localization-delocalization 
transition at K = 3/2. The derealization transition arises for much weaker attraction than 



in the one dimensional case [21] K p = 3. For the two chains problem the critical value 



of K can be attained for a Hubbard model p8| , |22[] whereas the one dimensional Hubbard 



model is always localized even for very negative U ||22j| . In addition the localization length 



2 1 

is increased : L 2c h. = {j)) z ~ 2Kp+ whereas in the one chain case L\ c h, = (jj) :i ~ Kp - At the 
opposite to what happens for the one dimensional case where the attractive localization 
length was smaller than the repulsive one pp|f2T[|, here the two lengths are the same, up to 



prefactors. Therefore, the enhancement of charge stiffness by repulsive interactions found in 



the one chain case [22] ought to be absent for 2 chains, or at least strongly reduced. Such 



an issue would need a more detailed study. The conductivity behaves as cr(T) ~ T 



2-2K 



p+ 



Clearly, these effects are due to the existence of a spin gap and to the freezing of interchain 
charge excitations |lTJ . As a consequence, it would be worth studying the localization effects 
in a three chains model (where there should be no spin-gap) to see if the derealization effect 
of attractive interaction does persist or if we fall back to the one chain case. 

We are grateful to H.J. Schulz for many useful discussions. 

Note added: After completion of this work, we learned about the work of Kawakami and 
Fujimoto |3"0fl . These authors considered the related, albeit different problem of disordered 
coupled Hubbard chains with a ferromagnetic Hund's exchange and no hopping. They also 
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find reduction of the localization effects in this system. 
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